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We discuss the inverse /3-decay of accelerated protons in the context of neutrino flavor superpo¬ 
sitions (mixings) in mass Eigenstates. The process p —>• n£ + vt is kinematically allowed because the 
accelerating field provides the rest energy difference between initial and final states. The rate of 
p —»• n conversions can be evaluated in either the laboratory frame (where the proton is accelerating) 
or the co-moving frame (where the proton is at rest and interacts with an effective thermal bath 
of £ and vi due to the Unruh effect). By explicit calculation, we show the rates in the two frames 
disagree when taking into account neutrino mixings, because the weak interaction couples to charge 
eigenstates whereas gravity couples to neutrino mass eigenstates pj. The contradiction could be 
resolved experimentally, potentially yielding new information on the origins of neutrino masses. 
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I. INTRODUCTION 

Neutrino mixings are an important piece of evidence 
for beyond standard model physics. The effect is phe¬ 
nomenologically described by the fact that neutrinos in¬ 
teract weakly with other standard model particles in fla¬ 
vor eigenstates \vf) (£ = e,/z, r) that are superpositions 
of mass eigenstates | vf) determined by a mixing matrix 
@, 1 : 

\ v t) =^2,Uti\vi) ( 1 ) 

i 

with Uu known as the Pontecorvo-Maki-Nakagawa- 
Sakata (PMNS) matrix. However, the fact that the mass 
basis and interaction basis cannot be simultaneously di¬ 
agonalized means the neutrino weak current and energy 
momentum tensor do not commute, and hence neutrino 
number and mass-energy are not commuting observables. 
This poses questions for their coupling to gravity M- 

To explore the neutrino-gravity coupling, we can per¬ 
form gedanken experiments involving neutrinos in accel¬ 
erated frames, using the fact that gravity is represented 
as invariance under general coordinate transformations, 
under which h for arbitrary 
infinitesimal vectors e M [§f. Specifically we can com¬ 
pare conventional Minkowski space quantum field the¬ 
ory calculations to curved-spacetime calculations by not¬ 
ing that transforming the Minkowski metric by = 
a -1 (ar, 0,0,1) yields the Rindler metric, representing 
an observer undergoing constant acceleration a. One 
might expect to obtain the same observables since the 
Minkowski and Rindler frames are related by a symme¬ 
try transformation of the theory, given that electroweak 
theories can be rewritten in a generally covariant way 
@, H3, H3, HI] • However, neutrino mixing violates condi¬ 
tions of the construction because, by breaking the mass 
superselection rule, neutrinos are not representations of 
the Lorentz group with a well-defined invariant P 2 . By 


considering in detail the results of calculations in differ¬ 
ent frames, we may be able to gain information or design 
experiments providing information on the physics of neu¬ 
trino mixing [l|. 

We consider the inverse /3-decay process p nv e e, be¬ 
cause it is potentially experimentally observable, but the 
effect is general and applies to any weak scattering occur¬ 
ring in an accelerated state. In the “laboratory” frame 
in Minkowski coordinates, the proton is accelerated by 
an external held and emits the electron and neutrino 
P4n|. The energy for the process is provided by the 
accelerating held, and the interaction is the electroweak 
vertex producing neutrinos in havor eigenstates. In the 
co-accelerating frame in Rindler coordinates, the proton 
is at rest and interacts with neutrinos and electrons in 
Rindler states [13, [Hj] , which display an effective thermal 
weight fl4| and are mass eigenstates, as we will show. 
However, mass eigenstates are related to havor eigen¬ 
states by Eq. [TJ Because the neutrino havor eigenstate is 
hxed by the electron l = e, we are forced in this frame to 
include the PMNS matrix in the amplitude. The PMNS 
matrix cannot be factored out, and the rates in the two 
frames differ. In this way, basic questions in quantizing 
neutrino helds Em are highlighted in the context of 
quantum held theory in curved spacetime IS in¬ 
consistency between coordinate frames has been ex¬ 
plicitly verihed in several cases: the Sokolov-Ternov ef¬ 
fect [19M21I ] . the emitted power in classical electromag¬ 
netic radiation [22|, E| and the rate of p —> nve + in the 
absence of neutrino mixing [ffl- Of these, only the 
Sokolov-Ternov effect has been seen in experiment [20j . 
though it may now be possible to look for the Unruh 
effect in the electromagnetic radiation of electrons ac¬ 
celerated by high intensity lasers [24], HU. Equivalence 
of observables is expected because the Unruh effect can 
be seen as a “Coriolis force” of quantum held theory in 
an accelerated state [2(| : an artifact of the coordinate 
system necessary to restore consistency, its derivation re- 
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quiring only Lorentz symmetry and quantum mechanics 

Our calculation implements the same technology for 
treating quantum processes in accelerated frames as pre¬ 
ceding work [9H 111] . We show that the limit of vanish¬ 
ing neutrino mass or trivial mixing Uu —> Sa recovers 
these results and the disagreement between frames arises 
only upon introducing the mixing. Although the mathe¬ 
matical origin of the disagreement is easily traced to the 
noncommutativity of weak and energy-momentum cur¬ 
rents, to understand the physics and use the experiment 
to learn something about neutrino mixing, we study the 
physics conditions required for the calculation to be valid. 

The equivalence of observables is shown to hold in the 
small acceleration limit, where the response of the sys¬ 
tem is linear and calculated from the classical state. We 
must be able to neglect both classical and quantum back- 
reaction on the accelerating field, such as radiation reac¬ 
tion in the electromagnetic case [23, [30] (beyond the lin¬ 
ear reg ime the effect of horizons needs to be considered 
[3 ll i3‘il |). This takes the form of a semiclassical approxi¬ 
mation, and consequently, the Hawking-Unruh effect has 
been compared to the Schwinger mechanism of particle 
production in classical electromagnetic fields [17j, [3J, |35|] ■ 
By detailing the approximations required to match the 
inertial frame result to the accelerated-frame result, we 
highlight several important differences, in particular how 
the details of the accelerating field are lost in the Unruh- 
type calculation. 

In the semiclassical limit gravity, described by the clas¬ 
sical metric perturbation field h M „, is required by Lorentz 
invariance to couple only to the conserved energy momen¬ 
tum tensor T [8|. As usual, we assume the neutrino is 
well-described as a free particle plus interactions, mean¬ 
ing that in the limit of no weak interactions its energy 
momentum tensor is that of three free fermion fields in 
their (ie mass) Eigenstate 


rpj_iV _ 


6L 


3 

y, | ViiD^^Vi + g^Vi(ilt) - m^Vi + ... 
i =1 


( 2 ) 

where D is the covariant derivative, comprising a spinor 
affine connection defined in [39|. Eq. [2] is diagonal in the 
mass basis, with mass eigenstates denoted by i = 1, 2, 3. 
At low energies -C Mw , the neutrino interactions are 
well-described by the Fermi effective theory. In all, the 
effective Lagrangian is 


3 

L = y ^(iy^ZF - nm)vi (3) 

i=l 

+ - me)tpe + y G F J Lh ^J^ t 

t £=e,n,r 


where J F is the left-handed current with nucleons h = 
p, n for proton and neutron. For leptons l = e, p, r, the 
weak-charge current is related to the neutrino mass eigen¬ 


states by 


Jli = y u k ( 4 ) 

i 

G f = 1.16 x 10 -5 GeV -2 is the Fermi constant. Its equa¬ 
tion of motion would show the metric field is sourced by 
the expectation value (T 7 *"), but here the kinetic term 
for the metric is neglected since we take the semiclassi¬ 
cal approximation and backreaction is negligible as long 
as the acceleration remains small. The lepton modes are 
determined including the classical field to all orders by 
solving the Dirac Hamiltonian, which is identical to the 
energy in a local rest frame, T 00 . The weak current is 
conserved, but, because of the PMNS matrix, it does not 
commute with T generating a potential source of ten¬ 
sion with gravitational physics even in the semiclassical 
limit0. 

In the following, we discuss the calculation in each 
frame, drawing attention to the approximations involved 
in the low acceleration limit and the necessity of the mass 
eigenstates and PMNS matrix in the Rindler frame. We 
shall perform the calculation in 1+1 spacetime dimen¬ 
sions for notational simplicity; the additional transverse 
directions do not affect the reasoning for needing the 
PMNS matrix in the Rindler frame. Then in Sec. ED 
we discuss experimental possibilities for resolving this 
paradox, and in Sec. lIVI we conclude with some broader 
lessons from this study. 


II. NEUTRINO MIXING IN DIFFERENT 
FRAMES 

A. Inertial frame small a expansion 

In the Minkowski frame (Fig. [TJ, the calculation 
amounts to electron and neutrino production by a classi¬ 
cal source [36|, because the Fermi theory current-current 
interaction is treated with a classical hadronic current 
J^e^Lh^ -A Jjj-^Lh n 0- The calculation is “semiclas¬ 
sical” in that quantized lepton fields are produced by a 
classical source. The time dependence of the source 
includes an oscillating phase gd" 1 **-" 1 ?) 4 contributing to 
the energy of the outgoing particles. To understand bet¬ 
ter what the process can reveal about the Unruh effect 
in presence of neutrino flavor mixings, we examine the 


1 In ordinary quantum mechanics, experiments such as Stern- 
Gerlach show that mutually non-commuting interactions do not 
give these problems, because different boundary conditions are 
generally related by rotations in Hilbert space. In quantum 
field theory, as we will see in the two subsequent sections, dif¬ 
ferent boundary conditions give unitarily inequivalent represen¬ 
tations. The tension between the gravitational and weak in¬ 
teraction terms, analogous to frustration in condensed matter 
systems, is therefore apparent only with a full quantum field 
theoretical treatment. 
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FIG. 1. Diagram for the process p —>• n£ + i'e in the inertial 
frame, with the wavy line corresponding to the classical po¬ 
tential, treated to all orders as indicated by the ellipsis. 


physics conditions necessary to reduce the hadronic cur¬ 
rent operator to a classical current, and this begins with 
an expansion in the small ratios of the acceleration to 
the other momentum scales a/mw,CL/m p ,a/ Am, where 
Am = m n — to p . 

To consider the proton motion prescribed, the acceler¬ 
ation must be small a -C m p , so that radiation reaction 
and backreaction on the external field is negligible. In 
fact, the proton’s compositeness gives a more stringent 
constraint since electromagnetic fields eB, eE ~ m 2 * af¬ 
fect its structure, so we require a <C m^/m p ~ 19 MeV. 
This implies also a -C mw so that the Fermi effective 
theory is applicable. The limitation to a/Am -C 1 will 
become clear later. The perturbative scale is a and hence 
the integration volume is d 4 x ~ a -4 , consistent with be¬ 
ing the physical scale of the saddlepoint in the tunneling 
potential for the leptons, which is localized (at the scale 
a -1 ) to the region around the turning point @ For the 
same reason, the acceleration should be constant over 
~ a -1 in order for the semiclassical approximation to 
apply for the leptons, which have be Broglie wavelengths 
of order Am. 

The p —> nve + operator from the lagrangian interac¬ 
tion in Eq. (0 

0 = ^^ n rPL%)(D elp P L e), Pl = ^^ (5) 

is a perturbative interaction, and the Fermi theory suf¬ 
fices since the momenta involved are all <C mw- To 
match the Rindler frame calculation which involves a 
correlation function restricted to a classical trajectory, 
we must reduce the proton to a point-like classical par¬ 
ticle. This differs from the semiclassical approximation 
in field theory, which entails solving the Dirac or Klein- 
Gordon equation in the presence of the classical accelerat¬ 
ing potential [37':, [Hj . To reduce to a classical trajectory, 
the anti-particle components of the proton wavefunction 
are integrated out first by using translational invariance 


2 In any case, at longer wavelengths the proton motion is per¬ 

turbed by the emission of radiation, which according to classical 

formulae has a wavelength of order a -1 . 


of the quasi-constant external held to focus on the re¬ 
gion around the turning point in the trajectory where 
the proton is nonrelativistic. Since a <C m n , the hadron 
is nonrelativistic for a time of order a^ 1 lnm„/m T > a -1 , 
found by solving for the time for the proton to gain lm p 
of energy. In this frame, the proton and neutron are 
heavy fermions, allowing us to extract the mass as the 
large part of the phase and integrate out the antiparticle 
components of the spinor 

E h P (v)e~ im P, h n(v')e imnt (6) 

V v' 

Corrections to the dynamics of the h v fields as well as 
antiparticle components are suppressed by l/m p . Also 
at leading order in 1 /m pi we can neglect the recoil, which 
means we drop the residual nucleon momenta and set v = 
v' . Considering the “coherence time” of the process to be 
St ~ 1 /mw, the velocity change during the interaction is 
Sv = aSt ~ a/mw , also a subleading correction. 

As the nucleons are nonrelativistic particles with neg¬ 
ligible recoil, the matrix element of the current satisfies 
the classical equation of motion in the external potential 
and therefore matches to the classical trajectory £ p (r) 
(essentially the Ehrenfest theorem) 

e -^mt {n \ K{v )rP L h P {v)\p) J^ cl {x) (7) 

with 

<d(*) = \ e ~ iAMt ^§) 5 ^ ~ &M) ( 8 ) 

The charge is included in the coupling constant Gf fac¬ 
tored out of the fermionic currents. Here, the factor 1/2 
comes from projecting an unpolarized current onto left- 
handed states. 

Evaluation of the matrix element proceeds straightfor¬ 
wardly, and the differential rate is 

dlFMink _ q 2 f ^ T ^ T ' e iAM(t-t')+i(ui e +u>l)(t-t')-i(k+k')-(^-i') 

d?kd 2 k' F J 

x e(k p kl + k'^kv - g pv k ■ k')u^ (t)u v {t') 

+ terms odd in k, k' (9) 

We have used dt/u 0 (r) = dr. We will discuss the phase 
space below, which comprises dependence on the final 
state lepton masses. Terms odd in k, k! vanish under the 
phase space integral and will be dropped in the following. 
The accelerated trajectory is parameterized as 

£ m (t) = a -1 (sinhar,0,0,coshar), u p = (10) 

dr 

The reason for the finite result for a process seemingly 
violating energy conservation is in this noninertial tra¬ 
jectory: The process is kinematically allowed due to the 
acceleration of the proton, because the external potential 
provides the rest energy difference between initial and fi¬ 
nal states. As long as the potential is “weak” a <C Am, 
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the process is exponentially suppressed, just as is spon¬ 
taneous pair production in quasi-constant electric fields 

min. 

The phase factor depends on the coordinates x^ with 
spatial coordinates replaced by C( T ) of the trajectory. 
Expanding around the origin where the hadron is in¬ 
stantaneously at rest, r = t — Ax ~ t(l — at/8) to 
leading order in at. Since the phase factor enforces 
energy-momentum conservation, t ~ Am~ l and we have 
t — t' = t — t 1 to leading (zeroth) order in a/Am. 

Making the substitution t — tf — > t — r', the integrand 
has no dependence on f = r + r', and one can define a 
rate per unit (proper) time. Defining T = J df we have 


1 dW M ink _ 1 dW 
T d 2 kd 2 k' T d 2 kd 2 k' 


Gl 


dcr e i2a 1 ( fe o+fco) sinhcra 


x (k 0 k ' 0 + k z k ’ z ) 


( 11 ) 


where the rescaled variables are 


k 0 = k° cosh to — k z sinh ra 
k z = k z cosh fa — /c 0 sinh fa 


We note that k and k are related by a Lorentz transfor¬ 
mation and, provided the phase space is Lorentz invari¬ 
ant, the rate WMink in Eq. [11] is also Lorentz invariant, 
defined in terms of the conversion reaction in its proper 
time. The transformation properties of the rate are ex¬ 
plicit in T~ l and thus WMink/? 1 is directly comparable to 
the rate obtained in the comoving frame. The rate mea¬ 
sured in the lab frame is obtained by convolving with the 
proton Lorentz factor along its lab-frame trajectory [§J 
as in the consideration of decays high energy particles. 

We discuss the phase space integration and comparison 
to the Rindler frame rate after obtaining the Rindler rate. 


B. Comoving frame calculation 


The Rindler space field operators are thus expanded 


X(x) = J ^ ( e %UV ~ b ™XuA u ) + -A u )) 

(13) 


in terms of the eigenmodes 

xw + = iv 1/2 


l 

0 




Xu - = * 7 ° 7 2 X^,+ 


V 0 


= K i^±± (men), N U1 1 = ^^ cosh— (14) 

a z CL CL 


The Rindler frequency takes all real values with no mass 
gap, and the mass appears in the wavefunction. It is 
the intrinsic dependence of the wavefunction on the mass 
that makes the sum over mass eigenstates not factoriz- 
able. For each field, the electron and each neutrino mass 
eigenstate, a corresponding complete set of solutions cov¬ 
ers the right spacelike wedge. 

Forming the matrix elements for the processes Eq. 
we take the Minkowski in and out states, bringing in the 
Bogoliubov coefficients relative to the Minkowski particle 
ij) + and antiparticle modes 


Xu = OtuA + Pu^P 

^ruj/ 2a 

(2 cosh(7rw/a)) 1 / 2 ’ 


Pu 


^—'Kuj/a _ 

= e 'a. 


(15) 

(16) 


This shows that in the Rindler frame, the p —> n transi¬ 
tion is matched to three distinguishable processes 


(?) 

pe — 

> nv e , 


Tu 


(II) 

p D e — > n e, 

Tu 

(17) 

(III) 

peu e - 

—> n 

TU 



The coordinates in the comoving frame are given by 
the Rindler metric, which in the right wedge 2 > \t\ has 
the form ds 2 = g llv dx ll dx L ' = u 2 dv 2 — du 2 . The space¬ 
like coordinate u is related to the Minkowski coordinates 
by u 2 = z 2 — t 2 (0 < u < oo) and determines the mag¬ 
nitude of the acceleration of the trajectory relative to 
the Minkowski space. We assign a continuous momen¬ 
tum eigenvalue u> conjugate to the time-like coordinate 
v, and the corresponding wavefunctions are determined 
by solving the Dirac Hamiltonian |39) 

uip = id v ip = Hif’ = 7°(y 3 u 1 ^ 2 d u u 1 ^ 2 — mt)ip (12) 

for each lepton t = e, z^. The 7 are the Rindler frame 
Dirac matrices, satisfying { 7 /i , 7 i '} = 2 . Note that 
the Hamiltonian operator is identically the g = v = 0 
element of the energy-momentum tensor operator. 


corresponding to the absorption by a static proton of, 
respectively, an electron, a neutrino and both an electron 
and a neutrino from the Unruh thermal bath (Fig. [2j). 

Additionally, the wavefunctions Eq. (fTH) do not cover 
all of Minkowski space, even when combined with the 
corresponding wavefunctions in the left, future and past 
wedge. A complete set of wavefunctions, in terms of 
which the Minkowski modes can be expanded, requires 
defining (5-function sources on the lightcones t ± z = 0 
fldj . For neutrinos, these sources are clearly in the mass 
basis and show that the wavefunctions in the calculations 
are required to be on-mass-shell solutions to Eq. (THZ1) at 
the boundary of the Rindler wedge. For this reason, we 
take the observed, asymptotic states in Rindler space to 
be mass eigenstates, and this requires introducing a fac¬ 
tor of Ui e to rotate to the electron field appearing in the 
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FIG. 2. Diagram for the process p —> n£ + i/e in the comoving 
frame, where the electron and neutrino legs are connected to 
thermal distributions. 


operator. Therefore, the squared matrix elements are 

(18a) 

\iM n \ 2 = jr(*.*oEi^i 2£ ^(*-* , )i^ e) i 2 i a - l) i 2 

(18b) 

(18c) 

where the hadronic and leptonic tensors are 
J^(x,x') = S^S^e- lAM(v - v,) S(u - a~ 1 )S(u' - a” 1 ) 


(19) 

G^(x,x') = ^2[x {vi) %P L X {e) ]( x ) \x (vi) lnP L X {e) }\x') 

( 20 ) 

in Rindler coordinates. Note that y 5 = 

ie» VK \ = 7 5 - 

We note that if the U*,U e i interference terms between 
mass Eigenstates were also added, inertial and comoving 
calculations would match, as this would be equivalent to 
assuming that a charge Eigenstate is the asymptotic state 
also in the accelerating frame. This would be consistent 
with the standard inertial-frame model of neutrino prop¬ 
agation, in which the different mass-states co-propagate. 
However including the interference terms would violate 
the Kubo-Martin-Schwinger (KMS) definition of a ther¬ 
mal state of a quantum system [T3, S3 by adding coher¬ 
ent, off-diagonal correlations in the density matrix. Con¬ 
sequently, the accelerated neutrino vacuum state would 
not be thermal, contradicting the essential characteris¬ 
tic of the Unruh effect and its formal derivation within 
quantum field theory [27j, H3 • It is possible that nature 
realizes this case, seeing that exact satisfaction of the 
KMS definition has not been proven for general interact¬ 
ing theories EM3 but such an alternative would have 
much more general implications than this work, given 
that the essence of the Unruh effect is the thermality of 
the state. 

Excluding the U*jU e i interference terms, we have fol¬ 
lowed the current best prescription for neutrino-emission 
and accelerated-frame calculations: weak interactions are 


perturbative, satisfying the condition to construct the 
usual wavefunctions for neutrinos in the frame of an 
accelerated observer (the proton) as free-particle solu¬ 
tions to the Dirac equation Eq. m- As the wavefunc¬ 
tions Eq. m have non-perturbative dependence on m„, 
boundary conditions require the asymptotic states in the 
accelerated frame are on mass-shell and contradict the in¬ 
ertial frame model of neutrino propagation, which does 
not require asymptotic neutrino states to be on mass- 
shell and is an experimentally-verified effective descrip¬ 
tion. The calculations in the inertial and accelerated 
frames are independently clearly defined, but, as we see 
now, imply incompatible results. This incompatibility 
can be seen as the price of maintaining the KMS condi¬ 
tion for a particle in a mass superposition. 

Summing the processes /,//,/// and integrating, the 
result has an analytic form similar to Q, with the addi¬ 
tional sum over neutrino mass eigenstates, 


w _ G 2 F a _ 2 
IU Rind - e 


i J-oo a a a 

( 21 ) 


where the integrand is 


,uj m P m ,,.. , ,m P . l0l 




Re 


ts- ( rUe \2 7v- t rn vj 


( 22 ) 

2 


with <*/ = ui — Am and K l ,(z ) are the modified Bessel 
functions of the second kind. The m Vi dependence of 
the integrand prevents factoring out the PMNS matrix 
and using )>T U*jU e i = 1. Without neutrino mixing, this 
result reduces to Eq. (9) of 0. 


C. Total rate comparison 

The simplest observable is the total number of 
electron-neutrino pairs detected, so that the Minkowski 
differential rate Eq. m must be integrated over the to¬ 
tal interaction time the electron and neutrino momenta 
in the final state. As usual, we ensure the 4-momentum 
integrals converge by putting the electron and neutrino 
on mass shell. Since the energy of the final state neutrino 
is E v ~ Am ^$> m„, it is a good approximation to the 
neutrino phase space to treat it as massless, 

WMi„k = f ^ 6(0 - m 2 e )5((k') 2 ) (23) 

The result converges to the m v —»• 0 limit given in [lT| 

WW = f duiH(—, ^,0) (24) 

87 t z J a a 

with H (x, y , z) the same function defined in Eq. E^l) . Ob¬ 
taining this form requires the integral transformations in 
M and the identity for the Meijer G-function 
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ofK v [x)K„{x) = ^GH 



i(l/ + /i + cr), 


These results carry through also for non-vanishing neu¬ 
trino mass, and the results of I, EH are reproduced. 

The inertial frame rate MA/iink Eq. (l24l) differs from 
the comoving rate Wrund Eq. ED since m v enters the 
wavefunctions in the comoving frame. 

The corrections to the integrand in Eq. (HHI) due to the 
771 j/ dependent terms appear to be small, but their exis¬ 
tence suffices to preclude factoring out JT £/Af7 e j. Since 
elementary particles live in different irreducible represen¬ 
tations no matter how small their masses, the identifica¬ 
tion of the created particle as a superposition of mass 
eigenstates survives to the infrared limit E 1 / 3 —> oo 
where the phase space is defined. The wavefunctions 
come with a factor E -1 / 2 (where E is the accessible con¬ 
figuration space volume) and the limit E —> oo must 
be taken with m„ finite for the approximation made in 
Eq. [23] to be valid. Taking m „ —> 0 first would remove 
the distinction between flavor and mass eigenstates and 
there would be no mixing, in contradiction to observa¬ 
tions but sufficient to make the inertial and comoving 
frame rates coincide. Therefore the massless limit does 
not commute with the infrared limit, and setting m v = 0 
in the Rindler frame Eq. [2T] would distort the infrared 
physics that identifies the asymptotic state as a mass 
eigenstate. This noncommutativity of limits is common 
to situations where infrared physics plays a crucial role 
(chiral symmetry breaking being an example Hi)- Its 
exact mathematical treatment most likely necessitates 
consideration of the issues described in the discussion 
section. 

On the other hand, one may worry that the neu¬ 
trino mass should have been included in the phase space 
Eq. (P251) . and doing so would lead to the inclusion of 
PMNS factors for the conversion to neutrino mass eigen¬ 
states. However, we recall that a neutrino is created in a 
flavor eigenstate, mixes over macroscopic distances and 
is detected in an experiment which again projects a flavor 
eigenstate. Since the neutrino energy is ~ Am, presently 
conceivable experiments would not be able to resolve the 
impact of the neutrino mass on the detection scattering 
event. This differs from particles in collider experiments, 
where detection occurs when the outgoing particle inter¬ 
acts with a macroscopic quasi-classical detector |44|j : the 
quantum evolution is effectively projected to a near-pure 
momentum eigenstate and reduces to a classical phase 
space integral in this limit [sj. 

As we pointed out in the introduction, this discrepancy 
is expected because of the ambiguity inherent in combin¬ 
ing non-mass Eigenstate fields with non inertial frames. 
In short, taking the m v —> 0 limit in the phase space 
in the Minkowski frame is a kinematical approximation 
separating corrections of order m^/E^ -C 1, whereas the 
TTiv —> 0 limit in the Rindler frame calculation makes an 


1 ^ i^,i 

2 cr, 2 (7 + 2 


- n + cr), i (—v + /r + cr), i(— v - n + a) 


(25) 


0(1) change in the infrared physics by removing mix¬ 
ing. Deforming to a nonflat metric deforms the mass 
eigenstate creation and annihilation operators. In the 
infrared limit, such operators define on-shell states, yet 
these are not selected by inertial detectors. Thus, one 
cannot expect descriptions in different frames to give the 
same scalar observables, as is generally true for interact¬ 
ing field theories. The resulting difference in the neutrino 
flavor distribution could in principle be observed by a 
flavor-sensitive neutrino detector ~ 1 oscillation length 
away from the accelerated-proton source, but we discuss 
potentially more achievable schemes in the section. 


III. EXPERIMENTAL OPPORTUNITIES 


This ambiguity warrants experimental investigation. 
A significant advantage to considering the inverse-beta- 
decay process is that it is reasonable to ask whether and 
how it could be experimentally verified. To test the above 
effect, we need only to accelerate protons and observe the 
rate of correlated neutron-positron emission. Although 
the physics is contained in neutrino fields, no neutrinos 
need to be detected. In contrast, weak processes with 
lower threshold, such as neutrino bremsstrahlung from 
an accelerated electron, require detecting the neutrino(s) 
in the final state. The authors know of no current plans 
for facilities where such an experiment is possible. Ac¬ 
celerated proton decay may be achievable with next gen¬ 
eration high intensity lasers. 

The challenges in observing this effect are (1) achieving 
sufficiently high a to overcome exponential suppression, 
and (2) maintaining that acceleration for sufficiently long 
that the quasi-constant approximation applies. The first 
is a matter of having a sufficiently strong and uniform 
field, something which is currently being developed in 
many contexts within pure and applied physics. The sec¬ 
ond may be partially offset by studying more general non¬ 
constant accelerated trajectories; however, in all cases the 
duration of the acceleration must be much larger than 
the “equilibration time” to the noninertial vacuum [45| . 
As a conservative estimate At < (am„) _1 ~ 6 fim=2 fs, 
though there is significant dependence on the acceleration 
profile. The field must also do work on the accelerated 
particles seeing that the final state rest energy is ~ AM 
greater than the initial state. 

Only electromagnetic fields can provide high accelera¬ 
tions for durations comparable to At. In the lab, these 
are provided by high intensity lasers. The current gen¬ 
eration of facilities, such as the Texas Petawa tt [4711 and 
the Extreme Light Infrastructure (ELI) [38j, |48(| create 
L ~ 20 — 40A ~ 20 — 40/im pulses and achieve fields of 
eE = (10“ 4 — 10 _3 )m 2 . A proton entering a field of this 
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strength is accelerated by a ~ 1CD 6 MeV, which means 
the rate oc e - 7rAm / a i s effectively zero. Next generation 
facilities are under discussion as part of a longer program 
to bring the QED critical field strength eE c = m 2 within 
reach of laboratory experiment. In a field of QED-critical 
strength, a proton’s acceleration is a ~ 10 -3 MeV, still a 
large suppression. If this process is to be experimentally 
realized, additional analytic study would be necessary, 
including likely nonperturbative sub-threshold enhance¬ 
ments, as have recently received much study in the con¬ 
text of spontaneous electron-positron pair creation, see 
for example f46| and related references. 

By comparison, other long wavelength fields are: (1) 
the mean electrostatic field inside an LHC bunch, esti¬ 
mated eE ~ (keV) 2 , using a cylinder with dimensions 
of the bunch (N = 10 11 , bunch radius 10~ 5 cm and 
length L ~ 10 cm); (2) mean fields inside fixed targets at 
hadron colliders eE ~ (100eV) 2 again with L ~ 10 cm; 
and (3) plasma wakcfield and crystal accelerator concepts 
eE < 10 _6 m 2 over L ~ pm. Larger L in the LHC bunch 
is not sufficient to overcome the relative e~ 10 suppres¬ 
sion. 

Quantum effects, in particular electron-positron pair 
production would be a large background. Although 
pair production is similarly exponentially suppressed, the 
lower threshold 2 m e < Am means it would occur at 
higher rate. For this reason, it would be necessary to 
solve the positron dynamics in the field to be able to cor¬ 
relate them with emitted neutrons. Scattering with elec¬ 
trons present in the system would also be a background, 
to be eliminated with ne + correlations and the known 
pe —> nv scattering cross-section. The theory work to 
understand these particle dynamics and backgrounds to 
greater precision is under way as part of the program 
to observe strong field QED effects in forthcoming laser 
experiments. 


IV. DISCUSSION 

We close with a discussion of what we expect to see 
if such an experiment is ever performed. The ambiguity 
between the inertial frame calculation and the comoving 
calculation Q, due to an interplay of mass and mixing 
(JJtj) terms in Eos 1^51 and [2T1 cannot be removed by a 
coordinate transformation. Either the comoving or the 
inertial calculation or both will give a wrong result, and 
we elaborate a few of the options revealed by making 
explicit the assumptions in the calculations. 

First, by writing the lagrangian Eq. (HU) , we consider 
the neutrino mass a tree-level operator, unmodified by 
the interaction producing the acceleration. This is con¬ 
sistent with the neutrino mass being either a “funda¬ 
mental” operator or generated by other interactions that 
have been integrated out, such as coupling to the Higgs 
or high scale beyond Standard Model physics. Measur¬ 
ing rates in agreement with the inertial frame calculation 
would provide evidence that the neutrino mass is an ef¬ 


fective operator of one of these types. Specifically, this 
treatment in the inertial frame is appropriate as long as 
the cutoff scale for the (effective) mass operator is a. 
Matching the comoving frame calculation to the observed 
rate in this case requires deeper understanding how the 
Higgs condensate transforms under the noninertial coor¬ 
dinate change or how effective operators with finite cutoff 
scales are represented in the curved space. 

It is also possible that neither calculation is an accurate 
prediction. Clearly, the model here is a crude approxi¬ 
mation to realistic experimental conditions using electro¬ 
magnetic fields to accelerate the proton: neutrons have 
zero charge and could not continue on the prescribed ac¬ 
celerated trajectory as the proton and the electron should 
be treated nonperturbatively in such a strong electro¬ 
magnetic field. In particular, conventional low energy 
(electromagnetic) initial and final state radiation can sig¬ 
nificantly modify the rate, needing to be resummed as 
is expected for high energy processes in high intensity 
lasers. 

An essential approximation is that the source of parti¬ 
cles is classical, though the radiated particles are treated 
as quantum, and we argued this approximation is valid 
when the acceleration is the smaller than any of the en¬ 
ergy scales associated with that current or the interac¬ 
tion. In this respect, the Unruh effect and specifically ac¬ 
celerated p —> nv e e + conversion are semiclassical calcula¬ 
tions and similar to spontaneous pair production in clas¬ 
sical electromagnetic fields. However, there is a crucial 
difference between the Unruh effect and the Schwinger 
effect: in the Unruh effect, the acceleration has been re¬ 
duced to a classical accelerated trajectory, and the de¬ 
tails of the accelerating field have been lost. All that one 
requires is a classical trajectory whose quantum fluctu¬ 
ations are minimized, so the action is real and close to 
the extremum. The Schwinger effect, on the other hand, 
necessitates an 0{K) correction to the action. 

This removal of the information about the acceleration 
is consistent with the original calculation in 0 which 
requires an essentially classical detector capable of pro¬ 
jecting single particle eigenstates. The nucleon here can 
be treated as a classical detector because distinguishing 
a neutron from a proton is “easy” for a classical macro¬ 
scopic detector, e.g. a particle tracker in a magnetic field. 

That said, both the Schwinger effect and the Unruh 
effect rely on the approximation that quantum fluctua¬ 
tions, controlled by H , are small. In the Schwinger effect, 
the occupancy number of the classical electromagnetic 
field is large 3> 1 implying that the field’s action is h, 
while in the Unruh effect it means a/8p h where 8p is 
the typical microscopic momentum exchange producing 
the acceleration. When the acceleration is produced by 
a field these two conditions are equivalent. On the other 
hand, too large a violates the assumption of no backre- 
action. 

The ambiguity highlighted here derives from breaking 
the mass superselection rule in the Lorentz group. Since 
the gravitational field must be sourced by the energy- 


momentum tensor, invariance of the h IJV T IJl/ coupling re¬ 
quires the particles states in the accelerated frame are 
mass eigenstates, which is the content of solving the Dirac 
equation in the Rindler coordinates Eq. (fl2l) . In the iner¬ 
tial frame, acceleration provides the energy required for 
the reaction but the interaction term always projects out 
flavor eigenstates. 

This ambiguity does not necessarily contradict the fact 
that the Unruh effect arises out of axiomatic field the¬ 
ory when written with general covariance [13, [13 ■ If it 
were possible to incorporate neutrino mixing into stan¬ 
dard quantum field theory with neutrino fields as well- 
defined representations of the Poincare group, then there 
would be a contradiction. However, the models of neu¬ 
trino mass and mixing in the literature all invoke inter¬ 
actions, for example with Higgs-like condensates or with 
high-scale particles integrated out. 

No axiomatic construction such as found in [13, H3 
is known for interacting theories, which are defined only 
perturbatively, so care needs to be taken in applying con¬ 
clusions from axiomatic field theories to theories with 
non-trivial vacua. It is well known that condensates 
break general covariance to zeroth order: in gravity, a 
Higgs-like mechanism is a common method to introduce 
consistently additional gravitational degrees of freedom 
[HI , or see [B(j for an example in the context of higher¬ 
dimensional extensions of the Standard Model. The full 
symmetry is restored only by constraining the higher- 
order interactions among new degrees of freedom on 
these backgrounds. Generally, low-energy effective ac¬ 
tions exhibiting neutrino mixing will not satisfy these 
constraints, and therefore we must determine the trans¬ 
formation of the condensates and interactions under gen¬ 
eral coordinate changes before we can say how the effec¬ 
tive interaction appears in the Rindler frame. 

Hence, the result presented here motivates systematic 
study of how effective interactions with high-scale par¬ 
ticles or condensates transform, or what conditions are 
necessary for the interactions to be consistent with gen¬ 
eral covariance. If neutrino masses originate in high- 
scale physics giving rise to low-energy effective interac¬ 
tions, the inertial frame tree-level calculation would be 
protected by the fact that the scale A responsible for 
the mixing (the Higgs condensate in the Weyl case, the 
non-renormalizeable operator in the seesaw case) is much 
larger than any momentum scale associated with the clas¬ 
sical field. The effective theory in which the inertial frame 
calculation was performed should therefore be correct to 
O (a/ A< 1). In the comoving frame, however, the tree 
level calculation would be based on an incorrect approx¬ 
imation since the neutrino would not be a point-like par¬ 
ticle but rather a composite between a “bare” degree of 
freedom and a zero-momentum condensate. In this case, 
the experimentally measured conversion would follow the 
inertial calculation, and an experimental realization of 
this paradox would teach us about the origin of the neu¬ 
trino mass. 

Another simple solution to the paradox is that the mix¬ 


ing matrix vanishes in the UV, just as for KaontJ^. In such 
a case, neutrinos with energy ~ Am could be mass de¬ 
generate. For consistency, also in the Rindler frame the 
neutrinos should lose the information of the mass eigen¬ 
state, and neutrino states with frequency^ Am appear 
as superpositions of the mass states. No PMNS matrix 
would be necessary, removing the contradiction. 

The proposed experimental program explores a physics 
domain orthogonal to that usually considered for beyond 
Standard Model physics: instead of one scattering event 
with large momentum transfers we consider a uniform 
acceleration, which means a large number of interactions 
with soft quanta. Since the effective potential giving rise 
to most, if not all mass terms is thought to be the Higgs 
mechanism involving a condensate of zero momentum 
quanta, it is reasonable to suppose new physics will show 
up in this regime provided the scale of the relevant term 
of the effective lagrangian is comparable to the accelera¬ 
tion, as it certainly is here. 

As a result of this orthogonality, the dependence on 
neutrino masses and mixing parameters of the rates cal¬ 
culated here can be very different from the usual ones. 
If the comoving frame calculation, Eq. [M] applies, it 
would open up an invaluable laboratory to study neu¬ 
trino masses and mixing angles directly: Unlike any other 
mechanism within known physics, the Rindler boundary 
conditions “prepare” the neutrino in a pure mass Eigen¬ 
state. The reaction rate depends on the absolute value 
of the phases within the mixing matrix and the abso¬ 
lute value of the masses, rather than mass differences. 
The experimental opportunities inherent in this realiza¬ 
tion are numerous enough that even enumerating them is 
beyond the scope of this work (For one, performing these 
experiments with both protons and antiprotons would 
give the CP-violating terms of the PMNS matrix). If the 
inertial frame rate Eq. [24] applies, the prospects of using 
the processes described here to study the PMNS matrix 
appear bleaker, since this matrix does not appear in Eq. 
[Ml However, precision studies of such a process might 
be able to shed light of the origin of the neutrino mass 
term, since loop corrections will certainly be sensitive to 
UV structure of the operator appearing as a mass in the 
IR. A detailed analysis of this is beyond the scope of this 
paper. 

In conclusion, we have argued that the conversion of 
accelerated protons into neutrons is a promising labora¬ 
tory to study the origin of neutrino masses and mixings. 
If the inertial frame calculation is correct, we do not have 
sensitivity to the neutrino mass absolute values, unless 


3 Kaons are mass-degenerate Eigenstates at the precision level the 
mass of an unstable particle can be conceivably measured l5ll . 
They oscillate because flavor does not commute with the weak 
isospin charge. Neutrinos oscillate because weak isospin charge 
does not commute with mass, as this work discusses. In the Kaon 
case, the mixing is due to an IR, operator coupling the Higgs 
condensate to Fermions, and it is possible a similar mechanism 
applies in the neutrino case. 
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our choice of the phase space is incorrect. Admittedly, 
technologically such experiments are still somewhat out 
of reach, but we hope that the ingenuity of the QED 
and intense laser community will make such experiments 
feasible in our lifetimes. 
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